Delay time and tunneling transient phenomena 
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Analytic solutions to the time-dependent Schrodinger equation for cutoff wave initial conditions 
are used to investigate the time evolution of the transmitted probability density for tunneling. For a 
broad range of values of the potential barrier opacity a, we find that the probability density exhibits 
two evolving structures. One refers to the propagation of a forerunner related to a time domain 
resonance [Phys. Rev. A 64, 0121907 (2001)], while the other consists of a semiclassical propagating 
wavefront. We find a regime where the forerunners are absent, corresponding to positive time delays, 
and show that this regime is characterized by opacities a < a c . The critical opacity a c is derived 
from the analytical expression for the delay time, that reflects a link between transient effects in 
tunneling and the delay time. 

PACS numbers: 03.65.Bz, 03.65.Ca, 73.40.Gk 



I. INTRODUCTION 

In recent times there has been relevant technological 
advances that have made possible to design and construct 
artificial quantum structures at nanometric scales jl], 
On the theoretical side the above achievements have 
stimulated work on the issue of time-dependent tunnel- 
ing. In particular, one finds a number of works that deal 
with the solution to the time-dependent Schrodinger's 
equation for cutoff wave initial states §, @ § §. One 
interesting feature of these approaches is that at asymp- 
totically long times the time-dependent solution goes 
into the well known stationary solution. This establishes 
a bridge between time-dependent and time-independent 
approaches that may be used to address some subtle 
questions, such as the controversial problem of the rel- 
evant time scales for tunneling 0. 

In a recent work we have used a time-dependent ana- 
lytic solution to the Schrodinger equation for an arbitrary 
potential^, to explore the tunneling dynamics for a rect- 
angular potential barrier ||. We found that the probabil- 
ity density exhibits a transient structure that we named 
time domain resonance, and obtained that it provides 
the largest probability of finding the tunneling particle 
at the potential barrier edge. Moreover we discussed the 
relevant time scales associated with the time domain res- 
onance as a function of the potential parameters and the 
incidence energy. 

The purpose of this work is to extend the above in- 
vestigation to study the time evolution of the probability 
density along the transmitted region of the potential. We 
found that for a large variation of potential parameters, 
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the probability density exhibits two evolving structures. 
One of them is a forerunner that corresponds to the time 
propagation of the time domain resonance, whereas the 
other structure consists of a propagating wavefront. We 
find that the forerunner vanishes at asymptotically long 
times and distances from the interaction region, whereas 
the propagating wavefront tends to the stationary solu- 
tion of the problem. The propagating wavefront exhibits 
a time delay with respect to the free case situation. Wc 
corroborate that this dynamical delay time is accurately 
described by the analytical expression obtained from the 
phase energy-derivative of the transmission amplitude. 
The analysis of this time scale as a function of the po- 
tential parameters yields positive and negative (time- 
advance) delay times. We have found that forerunners 
exist whenever the delay time is negative, thus establish- 
ing a deep connection between transient and asymptotic 
effects in tunneling. 

The paper is organized as follows: Section II provides 
the main expressions of the formalism that are relevant 
to calculate the probability density along the transmitted 
region. In Sec. Ill wc consider the time honored rectan- 
gular barrier potential model. Here we study through 
several subsections the time evolution of the transmitted 
probability density, and discuss the delay time. Finally, 
in Sec. IV we present the conclusions. 



II. FORMALISM 

The relevant expressions to calculate the time evolu- 
tion of the transmitted wave with the reflecting cutoff 
wave initial condition were considered in Ref. They 
follow from a general formalism developed by Garci'a- 
Calderon [^) for the solution of the time-dependent 
Schrodinger equation for tunneling through an arbitrary 
potential V(x) that vanishes outside a region < x < L. 
Our approach is a generalization of the free case problem 
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considered by Moshinsky || that led to the diffraction 
in time phenomenon. This transient effect has been re- 
cently verified experimentally^ and has stimulated fur- 
ther studies @. 

For the sake of completeness and to fix the notation 
we recall the relevant equations here. The cutoff wave 
initial condition corresponding to a reflecting wave may 
be written as, 



e- lkx , x<0 



ip(x,k;t = 0) 



(1) 



0, 



x > 0. 



The time-dependent solution ip(x,k;t) of Schrodinger's 
equation for the transmitted region, x > L reads, 

ip{x, k; t) = ip q (x, k; t) + xp r (x, k; t) (2) 

where ip q is given by, 

1> q (x,k;t) = T k M{y k ) - T_ fe M(y_ fe ) (3) 

and ip r by, 



ip r (x,k;t) = - T n M{y kri 



(4) 



In the above expressions the quantities T k and T_fc = Tj* 
refer respectively to transmission amplitudes and T n — 
2iku n (0)u n (L) exp(—ik n L) / (k 2 — is given in terms of 
the set of resonant states {u n (x)\ and the complex poles 
{fc„ = a n — ibn \ of the problem M. The functions M(y s ) 
are defined asH, 



M(y s ) = -^ mx ,2ht) w{iy s ), 



(5) 



where the w(iy s ) is the complex error function [JL2J with 
the argument y s given by 
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In Eqs. (||) and (^), s stands either for ±fc or k± n , and 
the index n refers to a given complex pole. Poles are 
located on the third and fourth quadrants of the complex 
/c-plane. The free case solution to the above problem for a 
cutoff plane wave was considered by Moshinsky ||. The 
solution for the free case with a reflecting initial condition 
is given by, 



ipa{x,k;t) = M{y k ) - M(p_ fc ). 



(7) 



Note that in the absence of a potential, i.e., T k = 1, 
the term tp q , given by Eq. (||), becomes identical to the 
free case solution ipo given above. We shall refer to 
ipq, that resembles the free contribution, as the quasi- 
monochromatic contribution and to the sum term given 
by Eq. (Q), namely i/v, as the resonant contribution. 



From the analysis given in Ref. pj one can see that the 
exact solution ijj = tp q + ip r , given by Eq. (||), satisfies 
the initial condition and that at asymptotic long times, 
ip r — * and ip q goes into the stationary solution. Hence 
at very long times tp becomes, 



ip(x, fc; t) = T fe c 



ikx —iEt I h 



(8) 



As pointed out in Ref. |J a cutoff wave initial state 
has, in addition to tunneling components, momentum 
components that go above the barrier height. One sees 
from Eq. (|^) , that the probability density exhibits and in- 
terplay between tunneling and over-the-barrier processes. 
However, as Eq. (||) indicates, at asymptotically long 
times, the transient effects vanish and one ends up with 
a stationary tunneling solution. 



III. THE MODEL 

As has been customary in studies involving tunneling 
times in one dimension we consider a model that has 
been used extensively in studies on time-dependent tun- 
neling, namely, the rectangular barrier potential, char- 
acterized by a height Vq in the region < x < L. To 
calculate the time-dependent solution ip(x, fc; t) given by 
Eq. (||), in addition to the barrier parameters Vo, L, 
and the corresponding incidence energy E = H 2 k 2 /2m, 
we need to determine the complex poles {k n } and the 
resonant states {u n (x)} of the system. Both the com- 
plex poles {k n } and the corresponding resonant eigen- 
functions {u n (x)}, can be calculated using a well estab- 
lished method, as discussed elsewhere. || ||. For the 
rectangular potential barrier the set of complex energies 
E n = h 2 k 2 /2m = e n — iT n /2 corresponds to the poles of 
the transmission amplitude of the problem^, and hence 
it may be used to describe the well known top-barrier 
transmission resonances appearing in that system. 



A. Dynamics of the transmitted probability density 

To exemplify the time evolution of the probability 
density in the transmitted region we consider a set of 
parameters typical of semiconductor artificial quantum 
structures §: Vo = 0.3 eV, L = 5.0 nm, E = 0.01 eV and 
m = 0.067 m e , with m e the electron mass. Our choice 
of parameters is the same as in Ref. ||, and it guaran- 
tees that most momentum components of the initial state 
tunnel through the potential. The different parameters 
may also be expressed in terms of the opacity a defined 
as, 



k L, 



(9) 



where fco = [SmVb] 1 / 2 /K, and the ratio u = Vq/E. In 
our case a = 3.63 and u — 30. For example, the regime 
of opaque barrier is reached for values of a > 5. In 
what follows we shall explore the time evolution of \ip\ 2 
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FIG. 1: Time evolution of the normalized probability den- 
sity \4>\ 2 (solid line) at the barrier edge x = L. The main 
contribution to the time domain resonance comes from the 
resonant term \ip r \ 2 (dotted line), the quasi-monochromatic 
contribution \ip q \ 2 (dashed line) oscillates with time in a sim- 
ilar fashion as the free solution \ipo\ 2 § (dashed-dotted line). 
The system parameters are: Vo = 0.3 eV, L = 5.0 nm, and 
E = 0.01 eV. See text. 

at several positions xq away from the interaction region. 
From Eq. (§) we can write, 

M 2 = M + \A\ 2 + i rq , (io) 

where I rq = 2Re(ip*4> r ) stands for the corresponding in- 
terference term. Figure [l] displays the time evolution of 
\ip | 2 (solid line) at the right edge of the potential bar- 
rier, x = L, as considered in Ref. This is the same 
example exhibited in Fig. 2 of Ref. |6| using a larger 
time scale. The sharp peak at very short times is mainly 
due to the resonant contribution \ip r \ 2 (dotted line). As 
discussed in Ref. Q, the maximum value of this time 
domain resonance, t p , provides a tunneling time scale 
representing the largest probability to find the particle 
at the barrier edge. In our example t p — 5.4 fs. The 
figure also shows the quasi-monochromatic contribution 
\ipq\ 2 (dashed-line) that rises and oscillates with time in 
a manner resembling the free case |"0o| 2 (dashed-dotted 
line). The interference contribution I rq is not shown, al- 
though clearly it is necessary to account for the complete 
solution. 

Along the transmitted region, x > L, the probability 
density becomes a propagating solution. We can see that 
the time domain resonance becomes a propagating struc- 
ture that we shall refer to as forerunner. Figure shows 
the case for xo = 50.0 nm. One sees that the amplitude of 
this transient structure (dotted line) is smaller than the 
quasi- monochromatic contribution (dashed line). Note 
also that the solution has separated itself into two well de- 
fined structures that propagate with different velocities. 
The forerunner propagates with a velocity given approxi- 
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FIG. 2: Time evolution of \i/j\ 2 (solid line) at the fixed po- 
sition xo = 50.0 nm. \ip r \ 2 (dotted line) and \i^q\ 2 (dashed 
line) account, respectively, by the forerunner and the quasi- 
monochromatic contributions. See text. 



mately by v r = hai/m, the velocity associated to the first 
top barrier resonance, whereas the quasi-monochromatic 
contribution does that, approximately by v k — Hk/m, 
the velocity associated to the incident particle. From a 
physical point of view we can understand the above sit- 
uation by noting that our initial state possesses momen- 
tum components in k-space above the barrier, which can 
be transmitted more effectively by the resonance window 
corresponding to the first top-barrier resonance. This is 
the origin of the fast tunneling response, given by the 
forerunner, mainly described by |Vv| • In our example, 
the main contribution to the resonant term, |Vv| 2 j comes 
from the first top-barrier resonant state. However, de- 
pending on the distance xq, one may need to sum up 
over many terms to account for the complete wave func- 
tion. The second type of response is given in a natural 
way by the maximum of the first peak of the probability 
density associated with the quasi-monochromatic contri- 
bution I'i/'ql 2 - This comes mostly from the momentum 
components centered about the momentum hk, which 
tunnel through the structure. 

At still much larger distances, xo — 1000.0 nm as 
shown in Fig. |3j the forerunner has disappeared almost 
completely, and the time evolution of the probability den- 
sity is dominated by |"0g| 2 - 

The behavior of the forerunner may be understood 
qualitatively by using the asymptotic properties j^] of 
the M(yk ±Tl ) functions in Eq. (||). By numerical inspec- 
tion we find that at a fixed position xq, the main fea- 
tures of the forerunner can be described using the one- 
term [n = 1) approximation to Eq. (Q), namely, \ipf \ 2 = 
\TiM (y^)] 2 ■ Since in the vicinity of the peak of the fore- 
runner the argument y^ of Mijj^) lies within —tt/2 < 
argy kl < tt/2, one obtains |?/;/| 2 = ^tt)" 1 \Ti/y kl \ 2 . 
This allows us to write a simple analytical expression for 
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FIG. 3: Time evolution of \i[>\ 2 at the fixed position xq — 
1000.0 nm. Notice that the forerunner has almost disappeared 
(an arrow indicates its position). The parameters are the 
same as in Fig. |l| 

the time evolution of this transient structure, namely, 

Vh? = ±_ (g/ro) 

|7i| 2 2tt [(*„ - fiaxt/m) 2 + {hht/m) 2 } ' 1 ' 

From the above equation we can see that the peak of 
the forerunner propagates with a velocity v r — Hai/m, 
as discussed earlier in the text. Figure ^ exhibits a plot 
of \ipf\ 2 (dashed line) as a function of time for the same 
parameters used in Fig. [I]. We observe good agreement 
with the exact calculation of \ifj\ 2 (solid line), given by 
Eq. (B). It is worthwhile to point out by inspection of 
Eq. (Ulh, that as time increases, the maximum of the 
forerunner, occurring at r = Xq/v t , diminishes at a rate 
proportional to Xq . Hence, for an increasing value of 
xq, the transient structure tends to a vanishing value. 

It is interesting to mention that in the case of opaque 
barriers, i.e., a > 5.0, the resonant contribution \ip r \ 2 
may be much larger than the quasi-monochromatic con- 
tribution I'f/'gl 2 . This occurs even at quite large distances 
from the interaction region. Figure |5j exhibits an exam- 
ple of this situation for L = 15.0 nm and a distance 
xq = 1 x 10 5 nm from the potential. Since the solution 
is normalized to \T^\ 2 and this quantity becomes very 
small for large L, one sees that \4>r\ 2 is several orders of 
magnitude larger than |?/> g | 2 , depicted in the inset to that 
figure. Clearly, as previously discussed, and exemplified 
in Fig. |^, at still much larger values of the distance Xo, 
the term \ifj q \ 2 shall eventually dominate over |^v| 2 . 

Let us now discuss another interesting behavior of the 
forerunner. In Fig. ^ we plot it as a function of time 
at a given distance xq > L, for different values of the 
barrier thickness, L = 5.0 nm (solid line), L = 4.5 nm 
(dashed line), L = 3.0 nm (dotted line), and L = 2.0 
nm (dashed-dotted line), for the same parameters used 
in Fig. pi We can see that the intensity of the transient 
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FIG. 4: Comparison between the formula for the forerunner, 
\ip/\ 2 (dashed line), and the exact solution, \ip\ 2 (solid line), 
as a function of time for a fixed value of the position xq = 50.0 
nm. The parameters are the same as in Fig. Ft]. 
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FIG. 5: The main graph shows the time evolution of \ip\ 2 
for the case of an opaque barrier of width L = 15.0 nm (a = 
10.88), at a fixed position xo = 1 x 10 nm. Notice that 
the forerunner, given essentially by the resonant contribution 
ji/vj 2 , overwhelms the quasi-monochromatic term described 
by IV 1 ?! 2 ) as depicted in the inset. See text. 



structure diminishes as L decreases. In fact, for the case 
of a barrier width L = 2.0 nm, we observe that the fore- 
runner disappears. However, as shown in the inset to 
Fig. @, what happens is that both the resonant contribu- 
tion, T^rl 2 ) and the interference term, I rq , in Eq. ( fio| ) for 
the probability density, are not only overwhelmed by the 
monochromatic contribution I'J'gl 2 but also almost cancel 
each other. 

Similarly Fig. [?] shows that the forerunner also disap- 
pears by diminishing the barrier height Vq , and again, as 
the inset shows, this occurs by the same reason as dis- 
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FIG. 6: Time evolution of the forerunner for a fixed value of 
the position xo = 50.0 nm, for different values of the barrier 
width L: (a) 5.0 nm (solid line), (b) 4.5 nm (dashed line), 
(c) 3.0 nm (dotted line), and (d) 2.0 nm (dashed-dotted line). 
Notice that the transient structure disappears as the barrier 
width diminishes. The parameters are the same as in Fig. ^ 
The inset exhibits a plot of the contributions to |^| 2 (dashed 
line) of case (d). Notice that the resonant contribution |^ r | 2 
(solid line) that gives rise to the forerunner, is almost canceled 
out entirely by the interference contribution I rq (dotted line) . 
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FIG. 7: Time evolution of the forerunner for a fixed value 
of the position xo = 50.0 nm, for three different values of the 
barrier potential Vb: (a) 0.3 eV (solid line), (b) 0.2 eV (dashed 
line), and (c) 0.1 eV (dotted line). Notice that the amplitude 
of the forerunner decreases as the barrier height of the po- 
tential diminishes. At the inset we plot the contributions 
to \fy\ 2 (dashed line) of case (c). Notice that the resonant 
contribution |* r | 2 (solid line) is almost canceled out by the 
interference contribution I rq (dotted line). The parameters 
are the same as in Fig. Ill 
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FIG. 8: Time advance of the solution \tp\ 2 (solid line) rel- 
ative to the free propagation case, | V^o | 2 (dashed line). The 
parameters are the same as in Fig. 1. 



cussed in the previous case. The above results hold also 
for the time domain resonance, i.e., at x = L. From the 
above analysis one could argue that the existence of the 
time domain resonance, and hence of the forerunners, 
depends basically on a particular combination of the pa- 
rameters Vb and L. In the next subsection we shall show 
that this is indeed the case. 



B. Delay time and forerunners 

An interesting result of the analysis of the previous 
subsection, depicted by Fig. ||, is that at very large dis- 
tances from the interaction region the time evolution of 
the probability density \ip\ 2 is essentially given by \^ q \ 2 
and exhibits a well defined wave front. As mentioned 
above, the wavefront propagates with approximately the 
classical velocity Vk = (hk/m), as follows by direct in- 
spection of the argument to the M— function, given by 
Eq. (^). A comparison of \ip\ 2 near the above wavefront 
with the corresponding free probability density \ipo\ 2 is 
exhibited in Fig. g. The parameters and the value of the 
position xo, are the same as in Fig. |[ Both solutions look 
very much alike. Note that its corresponding wavefronts 
are slightly displaced with respect to each other. In fact 
the maximum values of \ip\ 2 and iV'ol 2 ; for the parame- 
ters used in our example, exhibit a time difference that 
corresponds to a negative delay time (time-advance). 

The above considerations lead us to the notion of 
delay time as discussed by Bohm|l3|. He has argued 
that the main contribution to the transmitted proba- 
bility density comes from values in the neighborhood of 
space, for which the phase of the wave function changes 
slowly with energy. This yields the well known expres- 
sion for the delay iimejl3| as t^ = (d(f)/dk)/vk where <j> 
stands for the phase of the transmission amplitude, i.e., 
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FIG. 9: Plot of At (full squares) and the delay time t$ (solid 
line) as a function of the barrier width L, at fixed position 
xo — 1 x 10 J nm, with the same parameters as in Fig. 1. 
Here ko = 0.07258 and hence a varies from 7.25 x 10" 3 to 
13.065. See text. 



Tk = \Tk\ exp(i<j>) and Vk is the classical velocity as de- 
fined above. For the case of the rectangular potential 
barrier the delay time reads, 



m 



fc^sinh(2KL) - 2nLk 2 (k 2 - k 2 ) 



4k 2 n 2 + k* sinh 2 (kL) 



mL 
lik' 



. (12) 

We define At as the time difference of the maximum 
values of the curves for |V>| 2 and \tpo\ 2 obtained numeri- 
cally, namely, 



At = \^ max \ 2 - 



(13) 



Figure ^| displays a plot of At (full squares) and the delay 
time t,)) (solid line), as a function of the barrier width L 
for a large fixed value of the position xq . One sees that 
At reproduces exactly the behavior obtained from the 
analytical expression for t<$>. The above agreement of At 
with ttj,, does not hold when the distance xq is very close 
to the interaction region. There, the effect of the tran- 
sient structure cannot be ignored. Figure [l exhibits this 
situation for xq = L. The behavior of \ip 2 (solid line) 
is very different from that of the free contribution |f/>o| 2 
(dashed-dotted line). As discussed previously the time 
domain resonance peak comes from the resonant contri- 
bution and hence it is unrelated to the delay time t^. The 
splitting of the solution observed at larger distances and 
longer times has yet not occurred. 

Note also in Fig. ^ that for thin barriers there is a 
positive delay time, and as L increases, there is a tran- 
sition to a negative delay time. In what follows we shall 
demonstrate that such a transition occurs for a critical 
value of the opacity a c . In order to show the later, let us 
rewrite Eq. ( |l2| ) as a function of the parameters a and 
u = Vq/E, namely, 



[47 1 sinh 7 



cosh7 + 7 2 a 2 



-3] 



[7 2 a 



7 4 a 



74- 



sinh 2 ( 7 /2)] [4- 7 2 a- 2 ] 1/2 ' 
(14) 

where we have defined to = (mL/hko) and 7 = 2a(l — 
it -1 ) 1 / 2 . Thus from Eq. (|lj) the condition for the transi- 
tion from positive to negative delay times, i.e., t^ — 0, is 
simply given by the vanishing of the numerator, namely, 



47 sinh 7 — cosh 7 = 3 5- 



(15) 



For a particular value of the opacity a, we can deter- 
mine from the above equation the value of u at which 
the transition occurs. However, one finds by inspection 
of Eq. (p^[), that such a transition is not possible for 
small values of a. That is, there exists a critical value of 
the opacity a — a c such that for a < a c , the transition 
does not occur. This situation corresponds to impose the 
limit u — » 00 in the solution to Eq. (15). This implies 
that 7 — > 2a c, and hence Eq.(|l5|) becomes, 



cosh2cv r — 2a„ 1 sinh2a r = 1. 



(16) 



The numerical solution to Eq. ([16]) yields the critical 
opacity a c = 2.0653. Note that in addition to the po- 
tential parameters Vo and L, the opacity depends on the 
mass m of the incident particle. It turns out that this 
value of a c accounts for systems where the potential bar- 
rier is either too shallow or too thin. Therefore, in the 
regime a < cv c , only a positive delay time is observed. 

We have also found that a c plays an important role 
in the existence of time domain resonances and hence 
of forerunners. In fact, we find that for systems where 
a < a c , no time domain resonances nor forerunners are 
observed. This behavior can be seen in Figs. and [7]. 
For example, in the cases (c) and (d) of Fig. which 
correspond to barrier widths L = 3.0 nm and L = 2.0 nm, 
the parameter a is, respectively, 2.17 and 1.45. Clearly 
in case (d) the it forerunner has completely disappeared. 
In Fig. [?], the cases (b) and (c) corresponding to the 
potential heights Vq = 0.2 eV, and Vq — 0.1 eV, that 
refer, respectively, to the values of a, 2.96 and 2.095. In 
this case the disappearance of the forerunner occurs in 
the vicinity of the critical opacity a c . 



C. Comment on the phase-delay time 

Hartman has argued (ll| that the time it takes to a 
particle to traverse the classical forbidden region of a 
potential barrier, can be obtained from an analysis in- 
volving the delay time. He referred to this quantity as 
the transmission time th, though nowadays it is often 
called phase-delay time. It corresponds to the difference 
between the time at which a transmitted particle of mo- 
mentum hk would leave the rear of the barrier, x = L, 
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FIG. 10: Comparison of the transmission time th (solid line) 
with 5h (full squares) as a function of the barrier width L, 
measured at (a) the barrier edge x — L and (b) at a fixed 
position xo ~ 1 x 10 5 nm. The parameters are as in Fig. [l] 
Hence fco = 0.7258, and the opacity varies as mentioned in 
Fig. |. See text. 



and the time the same particle would arrive at the front 
of the barrier, x = 0. The transmission time th can be 
written as (see Eq. (13) of Ref. Q), 

th = U + t , (17) 

where t§ is given by Eq. (|l|) and t = {mL/hk) repre- 
sents the free time across a distance equal to the barrier 
width L. Note that to cancels out exactly the second term 
on the right-hand side of Eq. (|l7]) . The idea of consider- 
ing th as the relevant time scale for tunneling through a 
classically forbidden region has been criticized by argu- 
ing that there is no physical justification for relating in 
a causative sense the free evolving peak and the trans- 
mitted peak through a barrier [M Our analysis of the 
time evolution of the probability density supports this 
criticism. Indeed, as discussed in subsection A, along the 
transmitted region, the probability density may split into 
two structures evolving with different velocities. Hence 



it is not physically justified to choose a feature of one 
of them to compare it with the free evolving case. In 
particular, at the barrier edge x = L, for a > a c , the be- 
havior of probability density \ip(L,t)\ 2 is governed by the 
time domain resonance, that yields a completely different 
time scale|| than the phase-delay time. Moreover, even 
for a < a c , where there is no time domain resonance, our 
calculations do not support Hartman's transmission time 
th- This is illustrated in Fig. [To] (a), where we plot th 
as a function of the barrier width L (solid line) , and com- 
pare it with a plot of 5h = Ai+io (full squares) measured 
dynamically at the barrier edge x = L. Although both 
curves exhibit a similar qualitative behavior, the values of 
th and Sh, are quite different. On the other hand, Figure 
[To| (b) exhibits a plot of 5h (full squares) as a function 
of the barrier width L, measured at a distance xq very 
far away from the interaction region, i.e., xo = 1.0 x 10 5 
nm. This figure also shows a plot of th (solid line) and 
we see that they match quite well for all values of a. 
The lack of agreement between the plots of th and 5h in 
Fig. [l0| (a) follows because the time domain resonance, 
the quasi-monochromatic contribution and the interfer- 
ence term are very close together, as exemplified in Fig. 
ffl for L = 5.0 nm. On the other hand, at long distances 
the forerunner and the quasi-monochromatic contribu- 
tion are quite separated, though it may be shown that 
the interference term accounts for the delay timePq|. 

In the opaque barrier regime, a»l, the above times 
become independent of the barrier width, giving rise to 
the well known Hartman effect. Indeed at asymptotically 
large values of L, th goes as 2m/(hkK) as follows by 
inspection of Eq. ( fl7| ) pT| |. As can be seen, it is only 
at long distances from the interaction region, when Sh 
coincides with the dynamical time scale 5h- 

The above considerations, therefore, indicate that the 
transmission time, i.e., th given by Eq. (|l7|), does not 
represent the tunneling time of the particle through the 
classically forbidden region. The phase time delay t^, 
in the spirit of Wigner and Eisenbud|16 , represents an 



asymptotic effect of the potential on the tunneling parti- 
cle. 



IV. CONCLUSION 

Using an analytical solution to the time-dependent 
Schrodinger equation for cutoff semi-infinite initial waves, 
we have investigated the dynamics of the transmitted 
probability density for tunneling through a rectangular 
potential barrier. We have found two regimes, charac- 
terized by a critical opacity parameter a c , such that for 
values of a < a c there are no domain resonances and 
consequently no forerunners, whereas for a > a c , these 
transient structures may exist depending on the value 
of u = Vo/E. The above result follows from an unex- 
pected connection between the existence of these tran- 
sient structures and the delay time. This deserves to be 
further studied. An interesting feature of the formalism 
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used in this work is that it applies to arbitrary potential 
profiles of finite range. Hence, the existence of forerun- 
ners, for a given problem, would depend on the interplay 
among the different contributions to the probability den- 
sity given by Eq. ([To]). Our results suggest also that the 
study of transient effects cannot be ignored for a thor- 
ough understanding of the tunneling time problem. 
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